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Abstract: In two-Higgs-doublet models with exact Z2 symmetry, putting irih — 125 GeV at the align¬ 
ment limit, the following limits on the heavy scalar masses are obtained from the conditions of unitarity 
and stability of the scalar potential: vriH, rriA, < 1 TeV and 1/8 < tan/3 < 8. The constraints from 

6 —>■ 57 and neutral meson mass differences, when superimposed on the unitarity constraints, put a tighter 
lower limit on tan /3 depending on mfj+. It has also been shown that larger values of tan /3 can be allowed 
by introducing soft breaking term in the potential at the expense of a correlation between ruH and the soft 
breaking parameter. 


1 Introduction 


Extension of the Standard Model (SM) scalar sector is a common practice in constructing new physics models 
to address the shortcomings of the SM. Two-Higgs-doublet models (2HDMs) [1] are amongst the simplest of 
extensions that add only one extra SU (2) doublet to the SM scalar sector. The tree level value of the electroweak 
p-parameter remains unity for these types of extensions. In a general 2HDM, however, both the doublets ($1 and 
<i>2) can couple to each type of fermions. Consequently there will be two Yukawa matrices which, in general, are 
not diagonalizable simultaneously. This will introduce new flavor changing neutral currents (FCNC) mediated 
by neutral Higgses. It was shown by Glashow and Weinberg [2] and independently by Paschos [3] that Higgs 
mediated FCNC can be avoided altogether if fermions of a particular charge get their masses from the vacuum 
expectation value (vev) of a single scalar doublet. This prescription was realized by employing a Z2 symmetry 
under which one of the doublet is odd. Then there are four different possibilities for assigning Z2 parities to the 
fermions so that Glashow-Weinberg-Pascos theorem is satisfied. Following the usual convention, we shall always 
call <i>2 the doublet which couples to the up-type quarks. This leads to the following four types of 2HDMs: 

• Type I: all quarks and leptons couple to only one scalar doublet $2 ; 

• Type H: $2 couples to up-type quarks, while $1 couples to down-type quarks and charged leptons (minimal 
supersymmetry conforms to this category); 

• Type X or lepton specific: ^2 couples to all quarks, while $1 couples to all leptons; 

• Type Y or flipped: <i>2 couples to up-type quarks and leptons, while $1 couples to down-type quarks. 

Now that a light Higgs boson has been observed by the ATLAS and CMS Collaborations of the CERN Large 
Hadron Collider (LHC), the time is appropriate to revisit the constraints on 2HDM parameter space arising 
from the requirement of perturbative unitarity associated with the scattering amplitudes together with the 
complementary constraints coming from the flavor observables. Requirement of unitarity is a consequence of 
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probability conservation at the quantum level. Essentially, one takes tree level amplitudes of appropriate sets 
of scattering processes, and impose the unitarity condition that the s-wave scattering amplitude |ao| < 1. 
Equivalence theorem allows us to make the calculations easier by replacing the longitudinal gauge with the 
corresponding Goldstone bosons. In the context of the SM, Lee, Quigg and Thacker (LQT) had put an upper 
bound on the Higgs boson mass, trih < wlqt = (87r-\/2/3G_F)^/^ ~ 1 TeV [4]. The LQT-type analyses were 
later carried out to constrain the nonstandard parameter spaces for several extensions of the SM scalar sector. 
Detailed analyses of these constraints in the 2HDM context already exist in the literature [5-11]. In this paper, 
we study the consequences of adding two additional inputs to the existing analyses: (i) we use rrih — 125 GeV as 
input, which was an unknown parameter earlier, and also note that the fitted values of the couplings of the Higgs 
boson from the observed signal strengths into various fermion and gauge boson channels show close conformity 
to the alignment limit, and (ii) the restrictions coming from two crucial flavor observables, namely, branching 
ratio of 6 —>■ S 7 and the meson mass splittings have been superimposed on the unitarity constraints . We observe 
that for a large class of 2HDM with exact Zi symmetry, the unitarity and flavor constraints together allow a 
rather restricted zone for m}j+ and tan/3, thus imposing new constraints on these parameters. The limits on 
tan /3 are now independent of any other parameter, its upper limit coming from unitarity and the lower limit 
from flavor observables. We study these constraints in all four types of 2HDMs mentioned earlier. The limits 
are the strongest when the Z 2 symmetry is exact, while they get diluted by the soft symmetry breaking terms. 
We also comment on what happens when instead of Z 2 , softly broken U{1) symmetry is considered. 

This paper is organized as follows: in Section 2 we discuss the 2HDM scalar potential and the relevance of the 
alignment limit. Section 3 is divided into two parts. In the first part, we discuss the constraints arising from 
the requirements of unitarity and stability. Gonclusions obtained from this part are independent of the Yukawa 
structure of the model. In the second part, we revisit the Yukawa sector dependent constraints originating 
from flavor data and superimpose the result on the unitarity constraints. Finally, the important findings are 
summarized in Section 4. 


2 The scalar potential 


The general scalar potential of a 2HDM invariant under a Z 2 (4*2 —>■ —$ 2 ) symmetry can be written as [12] 


V = Ai 


A2(4>j4>2-^| + A 3 ( 4>[$i + $^$2 - 


Fit f 


Fit f 


+A 4 ((4>j;$i)(4'j4>2) - (4>];$2)(4-^4 >i)) + As (Rb 4*14*2 - + As (im 


4 * 14*2 


( 1 ) 


where, the bilinear term proportional to A 5 breaks the Z 2 symmetry softly. This soft breaking term, 
conventional parametrization, is written as —771^2(4*14*2 + 4>l4*i). The connection between 77732 A 5 is 
by [13] 

2 "*^5 

”^12 = ■ 


in its 
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( 2 ) 


Defining tan/3 = V 2 IV 1 to be the ratio of the two vevs, we also remember that it is the combination 
777^2/(sin/3 cos/3), not 777^2 itself, which controls the nonstandard masses [14[. In view of these facts, A 5 , rather 
than 77732 , constitutes a convenient parameter that can track down the effect of soft breaking. Note that, unlike 
the inert case, Eq. (1) implicitly assumes that the Z 2 symmetry is also broken spontaneously, i.e., both the 
doublets receive vevs. In this article, we shall only consider 2HDMs where the value of tan /3 is nonzero and 
finite. We have also assumed that all the potential parameters are real, i.e., GP symmetry is exact in the scalar 
potential. The last assumption allows us to define electrically neutral mass eigenstates which are also eigenstates 
of GP. Here there will be a total of five physical scalars: a pair of GP-even scalars [h and H with ran > rrih), 
one GP-odd scalar (H) and a pair of charged scalars (H^). For the transition into the mass basis, one needs to 
rotate the original fields in Eq. (1) by an angle /3 in the charged and the GP-odd sectors, whereas, the rotation 
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angle involved in the CP-even sector is denoted by a. The latter angle is defined through the relation 


^ o 2 (A3 -h V1V2 

tan 2a = --- , — -— . 

Xivf - X2VI + (A3 {vf - v^) 


( 3 ) 


Note that there were eight parameters to start with: vi, V 2 and 6 lambdas. One can trade vi and V 2 for 
V = \/vl + V 2 and tan /?. All the lambdas except A 5 may be traded for 4 physical Higgs masses and a. The 
relations between these two equivalent sets of parameters are given below : 
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(4b) 

(4c) 

(4d) 

(4e) 


Among these, v = 246 GeV is already known and if it is assumed that the lightest CP-even Higgs is what 
has been observed at the LHC, then mu = 125 GeV is also known. The rest of the parameters need to be 
constrained from theoretical as well as experimental considerations. 


The first simplification occurs if one keeps in mind that the experimental values of the Higgs signal strengths 
into different decay channels are increasingly leaning towards the corresponding SM predictions [15,16]. In the 
2HDM context, this implies the alignment condition [14] 


sin(/? — a) « 1 , 


(5) 


which means, h will have the exact same tree-level couplings with the vector bosons and fermions as in the SM. 
In view of the recent global fits for 2HDMs using the LHC Higgs data, Eq. (5) is a reasonable assumption [17-22]. 

Next, one has to ensure that there should not exist any direction in the field space along which the potential 
of Eq. (1) becomes infinitely negative, i.e., the potential is bounded from below. The necessary and sufficient 
conditions for this can be found to be [23,24] 


Ai -f A 3 > 0 , 

( 6 a) 

A 2 + A 3 > 0 , 

( 6 b) 

2 A 3 + A 4 -l- 2 \/(Ai -l- A 3 )(A 2 -f A 3 ) > 0, 

( 6 c) 

2 A,+ ^' + ''' l^= 2 ^'l+ 2 v'(A, + A,)(A, + A,)> 0 . 

( 6 d) 


To obtain the constraints from tree-unitarity, we construct an S'-matrix using different two-body states to label 
its different rows and columns. The £ = 0 partial wave amplitudes for different 2 —>■ 2 scattering processes 
constitute the elements of this S'-matrix. The explicit expressions for the eigenvalues of this matrix are listed 
below [5-8] : 
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Figure 1: Allowed region from unitarity and stability for exact Z 2 symmetry. In the right panel, the dark (purple) 
shaded region is excluded from direct search [25]. 
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= (Ai + A2 + 2A3) ± ^(Ai - A2)2 + ^ (A5 - Ae)' , 
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bb 
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be 
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( 71 ) 


The requirement of tree-unitarity then restricts each of the above eigenvalues as 

|af I, |6,| < 167r. (8) 

Now it is the time to investigate how the constraints from unitarity and stability restrict the parameter space 
in the alignment limit defined by Eq. (5) . 


3 Constraints on the 2HDM parameter space 


To begin with, we consider the case when the Z 2 symmetry is exact in the scalar potential i.e., A 5 = 0. For this, 
we have generated millions of random points in the {tan/3, mn, rriA, rnu+'\ space. The individual parameters 
have been varied in the following range: 

tan/3 G [0.1,100], mn G [125,2000], tua G [0,2000], G [0,2000]. (9) 

Those points which successfully negotiate the unitarity and stability bounds have been plotted in Fig. 1. Some 
noteworthy features are listed below : 

• From the left panel, one can read the limit on tan/3, 1/8 < tan/3 < 8 . 

• Strong correlation exists between the upper limit of mu and tan /3. 

• Limits on the masses are, rriA, rnff+ < 1 TeV. 
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Figure 2: Relaxation of the unitarity and stability constraints on tan ft in the presence of soft breaking. In the 
right panel, the vertical width of the of the tail in the region where tan (3 is much away from unity is caused by 
the variation 0 /A 5 in the range [—15,15]. 



Figure 3: Effect of soft breaking on the constraints on the nonstandard masses. The (light blue) shaded region in 
the rightmost panel represents the combined allowed region from direct search and the diphoton signal strength 
at 95% C.L. 


The reason for the above features can be traced back to the eigenvalues of Eq. (7). First two constraints for 
boundedness in Eq. ( 6 ) can be combined into 


Ai + A 2 + 2 A 3 > 0 . ( 10 ) 

This, then together with the condition |a]’^| < IGtt, implies 

0 < Ai + A2 + 2A3 < —^ , 

^ 0 < (tan^ P + cot^ /3) + 2m^ < 

where the last expression is obtained from the previous one by using Eq. (4) in the alignment limit. Keeping 
in mind that mn > 125 GeV, this will put a limit on tan/3 (as well as cot/3) when A 5 = 0. Since the minimum 
value of (tan^/3 + cot^/3) is 2 when tan/3 = 1, the maximum possible value of mn occurs at tan/3 = 1. In 
summary, the inequality (12) explains the tan/3 dependent bound on mn as depicted in the left panel of Fig. 1. 


( 11 ) 

( 12 ) 


Eq. (12) also implies that the restriction on tan/3 will be lifted for 1 / 2 A 5 U^ > (m|j)min = (125 GeV)^. Once 
this condition is satisfied, m'jj will have the chance to saturate to 1 / 2 A 5 U^ making the difference between them 
to vanish in Eq. (12). In fact, to a very good approximation, one can use 

mjj « I/ 2 A 5 U 2 (13) 


for tan /3 > 5. 
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To understand the restrictions on tjia cind Tn//+, we use the triangle inequality to note the following : 

1 

| 6 i - 63 I = 3|A6 - A 5 I < 327r, ^ \m\-< ——, (14a) 

1 1 ^TTl)^ 

|66 - = 3 |A 4 - A 5 I < 327r, ^ ^—■ (14b) 

Because of Eqs. (14a) and (14b) one expects to put limits on niA and m[{+ respectively, when A 5 = 0. Addi¬ 
tionally, note that due to the inequality 


107ru ^ 

|5i - bel = 3|A6 - A 4 I < 327r, ^ |TO^-m^+|< — - —, (15) 

it is expected that the splitting between rriA and TO 44 + will be always restricted in a 2HDM. It is also interesting 
to note that the conclusions obtained from Eqs. (14a), (14b) and (15) do not depend on the imposition of the 
alignment condition. 

Next we shall investigate the implications of the soft breaking parameter on these constraints. We have varied 
As in the range [—15,15] for this purpose. From Eq. (12) one can observe that the space for tan/I is squeezed 
further if A 5 < 0 but the bound is relaxed if A 5 > 0. This feature emerges from the left panel of Fig. 2. One 
can also see from Eq. (12) that must follow l/ 2 A 5 n^ if tan/3 moderately deviates from unity. This feature 
is reflected by the horizontal tail in the right panel of Fig. 2 on both sides of the peak. The vertical width of 
the tail is caused by the variation of A 5 in the range [—15,15]. On the other hand, from Eqs. (14a), (14b) and 
( 12 ), it should be noted that the upper bounds on the nonstandard scalar masses will be relaxed for A 5 > 0 but 
will get tighter for A 5 < 0. Fig. 3 reflects these features where one can see that this dependence is rather weak. 
It is worth mentioning at this point that if one uses a softly broken U{1) symmetry instead of the usual Z 2 
symmetry, the soft breaking parameter gets related to the pseudoscalar mass as m\ = l/ 2 A 5 n^. Consequently, 
the correlation between and A 5 in the leftmost panel of Fig. 3 transforms into the degeneracy between 
niH and niA- Detailed analysis of the scalar sector, for the softly broken 17(1) scenario, has been carried out 
in [26,27]. 

It is also important to note that the production as well as the tree-level decay widths of h remain unaltered from 
the corresponding SM expectations due to the imposition of alignment limit of Eq. (5) . But the loop induced 
decay modes of h, such as h ^ and h —>■ Z 7 , will pick up additional contributions due the presence of the 
charged scalar in loops. For example, the diphoton signal strength (/i-y-,-), in general, depends on both A 5 and 
771/4+ [13]. The current measurement by CMS gives = 1.14/(g’23 [28], whereas ATLAS measures to be 
1.17 ± 0.27 [29]. In addition to this, the direct search limit of TO//+ > 80 GeV [25] should also be taken into 
account. Considering all of these experimental constraints, the allowed region at 95% C.L. has been shaded (in 
light blue) in the rightmost panel of Fig. 3. Only those points that lie within the shaded region survive both 
the theoretical and experimental constraints. 


3.1 Yukawa sector and flavor constraints 


Now we shall concentrate on the constraints on the charged scalar mass (m/f+), imposed by the measured values 
of 5 —>■ 57 branching ratio [30] and neutral meson mass differences (AM) [31]. Since we are concerned with the 
quark sector only, the constraints will be the same for Type I and Type X models. The same is true for Type II 
and Type Y models. In the following, we spell out the relevant parts of the charged scalar Yukawa interaction: 


^(I or X) 

^(II or Y) 
“^ 77 + 


^/2H+ 

V 

y/2H+ 

V 


cot fd^UR {DuV) dL - UL (VDd) di/j -k h.c. , 

I cot/3 ur {DuV) di -k tan/3uL {VDj) d/jj -k h.c. 


(16a) 

(16b) 
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Figure 4: Constraints on tan /? and the charged Higgs mass from unitarity and flavor physics. The left panel 
corresponds to Type I and X and the right panel to Type II and Y scenarios. The lower horizontal dark (purple) 
strip in both the panels corresponds to the direct search limit of 80 GeV [25]. The ligther shades represent 
allowed regions from individual flavor observables. The scattered points are allowed from unitarity and stability 
for 2HDMs with exact Z 2 symmetry. 


where, V is the CKM matrix and D^.d are diagonal mass matrices in the up- and down-quark sectors respectively. 
In writing Eq. (16), we have suppressed the flavor indices. Thus u and d should be interpreted as three element 
column matrices. 

For the process b sj, the major new physics contributions come from charged scalar loops. We have added 
the new physics contribution to the SM one at the amplitude level and therefore have taken the interference 
into account. The branching ratio is then compared with the experimental value, (3.55 ± 0.26) x 10“^ [30], to 
obtain the allowed region at 95% C.L. in Fig. 4. As can be seen from Eq. (16b), for Type II and Y models, in 
the charged Higgs Yukawa interaction, the up-type Yukawa coupling is multiplied by cot j5 while the down-type 
Yukawa is multiplied by tan /3. Their product is responsible for setting tan /3-independent limit^ 'raH+ > 320 
GeV for tan/3 > 1 [33-35]. This feature has been depicted in the right panel of Fig. 4. In Type I and X models, 
on the other hand, each of these couplings picks up a cot (3 factor. This is why there is essentially no bound on 
mfj+ for tan > 1 in these models [34] . This character of Type I and X models emerges from the left panel of 
Fig. 4. 

The dominant new physics contributions to neutral meson mass differences come from the charged scalar box 
diagrams. Note that, the constraint arising from A/VL in the TOj:/+-tan/3 plane, is slightly stronger than that 
from the precision measurement of the Z ^ bb branching ratio [35] . In Fig. 4, allowed regions have been shaded 
assuming that the new physics contributions saturate the experimental values of AM [31] . Since the amplitudes 
for the new box diagrams receive prevailing contributions from the up-type quark masses which, for all four 
variants of 2HDMs, comes with a cot /3 prefactor, the overall charged scalar contribution to the amplitude goes 
as cot^ (3 due to the presence of four charged scalar vertices in the box diagram. Not surprisingly, AM offers a 
stronger constraint than b ^ sj for tan (3 < 1 because, in this region, the new physics amplitude for the latter 
goes as cot^ /3. 

Things become more interesting when the above flavor constraints are superimposed on top of the constraints 
from unitarity and stability. Most stringent constraints are obtained when Z 2 symmetry is exact in the scalar 
potential, i.e., A 5 = 0. In Fig. 4, the scattered points span the region allowed by the combined constraints 
of unitarity and stability for the case of exact Z 2 symmetry. Only those points which lie within the common 
shaded region survive when all the constraints are imposed. For 2HDMs of all four types, one can read the 

^Recently this bound has been changed to m^+ > 480 GeV using the NNLO result along with the updated experimental 
value [32]. This means that the b —>■ sy boundary on the right panel of Fig. 4 will be shifted upwards slightly. But this will hardly 
modify the lower bound on tan/3 which mainly comes from AM. 
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bound on tan j3 as 


0.5<tan/3<8. (17) 

However, in order to allow a lighter charged scalar in the ballpark of 400 GeV or below, one must require 
1 < tan/3 < 8. It should be remembered that, the lower bound on tan/3 mainly comes from the flavor data, 
whereas the upper limit, for the case of exact Z 2 symmetry, is dictated by unitarity and stability. In the presence 
of a soft breaking parameter, the upper bound will be lifted allowing tan /3 to take much larger values at the 
expense of a strong correlation between the soft breaking parameter and mu as depicted by Eq. (13). 


4 Conclusions 


In this paper, we have revisited the constraints from tree-unitarity and stability in the context of 2HDMs. The 
observed scalar at LHC has been identified with the lightest CP-even scalar of the model. The alignment limit 
has been imposed in view of the conformity of the LHC Higgs data with the SM predictions. These are the new 
informations that became available only after the Higgs discovery. If the Z 2 symmetry is exact in the potential, 
it is found that all the nonstandard masses are restricted below 1 TeV from unitarity with the upper limit on 
mn being highly correlated to tan/3. The value of tan/3 is also confined within the range 1/8 < tan/3 < 8 from 
unitarity and stability. The constraints from flavor data severely restrict the region with tan/3 < 1. Therefore, 
for an exact Z 2 symmetry, tan j3 is bounded within a very narrow range of 1 < tan f3 < 8 when a light charged 
scalar with mass around 400 GeV is looked for. 

In the presence of an appropriate soft breaking parameter the upper bound on tan /3 will be diluted. However, 
for large values of tan /3, the unitarity and stability conditions will render a strong correlation between the soft 
breaking parameter and mn as appears in Eq. (13). It is also worth noting that the value of can play a 
crucial role in the presence of soft breaking. For example, if Hjy is measured to be consistent with the SM 
expectation with 5% accuracy then one can conclude [13] for any value of tan/3. Thus, for 

large values of tan/3 one may expect ~ l/2A5n^. In this limit, the heavier nonstandard scalars 

truly decouple from the low energy observables. Thus, as has been emphasized in [13], proper decoupling of the 
nonstandard scalars necessitates the presence of a soft breaking term in the scalar potential. 

To sum up, when flavor constraints are superimposed on the constraints from unitarity and stability, the value 
of tan/3 is restricted within a very narrow range of 1 < tan/3 < 8 for 2HDMs with exact Z 2 symmetry. Larger 
values of tan /3 can be allowed by introducing suitable soft breaking parameter in the scalar potential, but the 
theoretical and experimental constraints impose certain correlations between nonstandard masses and the soft 
breaking parameter. This makes the theory much more predictive in the large tan/3 region. 
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